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We investigate classical-quantum correspondence for 
kicked Harper model for extremely small values of the Planck 
constant h. In the asymmetric case a pure quantum state 
shows clear signature of classical diffusive as well as super dif- 
fusive transitions asymptotically independent of h. However, 
for the symmetric case, the h independent behavior occurs 
only for renormalized parameter K — K/{2h) with intriguing 
features such as a sharp transition from integrable to non- 
integrable transport ai K = ■7'"/2, a series of transitions at 
multiples of n and periodicity of the transmission probability. 
These results add new puzzles to the frontiers of quantum 
chaos. 

PACS numbers: 72.15.Rn-h72.15-v 



Localized transport of quantum system in the regime 
where the corresponding classical system exhibits de- 
terministic diffusive behavior is one of the most sur- 
prising aspects of non-integrable Hamiltonian systems 
[n|. However, the correspondence principle requires some 
signatures of various classical transitions such as the 
breakup of KAM tori leading to diffusive transport and 
the emergence of accelerator modes (AM) resulting in 
super-diffusive anomalous transport. In this paper, we 
describe quantum signatures of various classical transi- 
tions in transport characteristics emphasizing the cross- 
over effects from large to small values of effective Planck 
constant h. As h —^ 0, the quantum system exhibiting 
localization in one of the phase-space directions is found 
to feel the effects of all classical transitions. However, 
in the absence of localized transport, the quantum sys- 
tem exhibits many surprising features and appears to be 
insensitive to the classical dynamics. 

The problem of establishing classical-quantum corre- 
spondence in quasi-periodic extended systems has proven 
to be difficult due to numerical limitations in approach- 
ing ?i "> 0. All previous studies (see e.g. p) addressing 
this question have been limited to ?i ~ 1. Here we use 
recently developed renormalization group (RG) approach 
to study quantum transport for extremely small val- 
ues of h upto ~ lO"**. 

The kicked Harper model |Q,|| has emerged as an im- 
portant model in quantum chaos literature. The system 
is given by doubly periodic time-dependent Hamiltonian 



Here q,p is a canonically conjugate pair of variables, usu- 
ally considered on a cylinder p G (— oo, oo), q G [0, 27r). 

The classical dynamics of the kicked Harper is deter- 
mined by two parameters K and L. In the asymmetric 
case {K ^ L), the phase space, for small values of pa- 
rameters, is stratified with KAM tori which inhibit the 
transport in global scale. For L > K, or K > L, these 
tori barriers limit the transport along p, or x, directions, 
respectively. As we show below, the 2-d parameter space 
exhibits an intricately mixed non-diffusive (KAM regime) 
and diffusive regions corresponding to global stochastic- 
ity. In contrast, in the symmetric case K — L, there 
are no KAM barriers to global transport but a separa- 
trix and it is the breakup of the separatrix that results 
in global diffusion. For large values of the parameters, 
both the symmetric and the asymmetric model exhibits 
mostly diffusive behavior with the exception of narrow 
windows in parameter space where the AMs give rise to 
super-diffusive transport, [g 

The quantized system that is periodically kicked is de- 
scribed by the quasi-energy states of the one step time 
evolution operator, introducing an additional parameter 
h into the problem. However, one hopes to recover h 
independent behavior (as 7i ^ 0) in order to establish 
quantum signatures of classical behavior. It is a well es- 
tablished fact that the RG approach provides the most 
effective tool in distinguishing ballistic, diffusive and lo- 
calized transport. Here we use recently developed O] 
dimer decimation approach to study transport charac- 
teristics of the quasi-energy states in the small h limit. 
The RG method is applied |^ to the momentum lattice 
(Pm = hm) representation of the kicked model 
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where the coefficients B™ are 

B™ = Jr (K) sin[L cos(to?i) - 7rr/2 - u/2] . (3) 

We introduce renormalized parameters as K = K/(2K), 
L = L/{2h). Tight-binding model (TBM) (|) effectively 
contributes only few terms as Bessel's function exhibit 
fast decay when \r\ > \K\. Therefore, the TBM describes 
a lattice model with a finite range of interaction denoted 
as b {b K^ K). In the limit of small K,L,uj, TBM re- 
duces to the Harper equation with e = huj j|]. We will 
choose h to be an irrational number with a golden tail: 
h/{27T) = l/(7ih + cr) where a = (\/5 - l)/2) is fixed 



and h is varied by varying the integer n\^. This corre- 
sponds to studying system sizes Nn,n — 1,2... deter- 
mined from the Fibonacci equation iVn+i — ^n + -^n-i, 
with A^o = IjA^i = »T-h, and corresponding to n-th suc- 
cessive rational approximant of the irrational number g. 
The RG methods can be used to study system sizes upto 
10^ which allows very large nh and hence facilitates 
studying kicked model for extremely small values of fi. 

The transport characteristics of the quasi-energy states 
are studied by computing the transmission probability T 
on the momentum lattice. This is achieved in two steps: 
first we decimate the lattice and then solve the scattering 
problem on the renormalized lattice P| . Renormalization 
scheme makes the solution of the scattering problem for 
large lattices of size N very efficient, as the dimer dec- 
imation reduces the size of the scattering region. For a 
fixed ?i, we compute the transmission probability T{N) 
for various sizes N of the momentum lattice correspond- 
ing to a rational approximant of a with denominator N. 
The scaling exponent (3 — limjv^oc log T(iV)/ log iV dis- 
tinguishes extended, localized and critical states respec- 
tively, described by fi{N) — > 0, ^ — oo, and oscillatory 
function f3{Nn) of n B. For the exponential localization, 
the quantity ^ = -[limN^oail/N) \ogT{N)]~^ has been 
found to be closely related to the localization length of 
the quasienery eigenstate to. 
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nent /3c\ = lim^^oo log((p(i) - p{Q))'^)/\ogt, are signaled 
by /3ci changing from to ~ 1. In the quantum model, our 
detailed analysis for various values of H and L > K, con- 
firms the previously held view |^ that the quantum sys- 
tem remains localized in the classically diffusive regime. 
However, the classical transitions corresponding to diffu- 
sive transport manifest in huge enhancement of localiza- 
tion length. Results for an individual pure state uj — 
are shown in Fig. 1. As demonstrated in the figure, small 
h is crucial to see signatures of all classical transitions. 
We would like to point out that our results are consis- 
tent with the relation ^ — ^D/h? |0]. However, near the 
peaks, (narrow windows in parameter space correspond- 
ing to the onset to classical transitions), the quantum 
transmission probability T{N) ceases to look like a sim- 
ple exponential ~ exp(— A^/^) thus making quantitative 
comparison difficult. 

An interesting aspect of two-parameter Harper kicked 
model is that the boundary between the KAM and diffu- 
sive phases appears to be fractal as seen in Fig. 2. This 
behavior is reminiscent of the kicked rotor problem where 
the kicking potential consists of two-harmonics |l3] . Al- 
though somewhat smeared, the quantum model exhibits 
similar behavior: the boundary describes the transition 
to the enhancement of localization length ^. It is re- 
markable that the quantum system feels the presence of 
all classical transitions and the fact that unlike kicked 
rotor, there is a whole hierarchy of transitions in Harper 
model, makes this model an important system in quan- 
tum chaos studies. 
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FIG. 1. For fixed K — 0.4, the figure describes a series of 
breakup of KAM barriers resulting in transitions to difi'usive 
transport /3ci = 1 (top). Middle and the bottom figures show 
the corresponding plots of quantum localization length ^ for 
rih = 200 and rih — 32. The renormalization is carried out 
for system sizes increasing until the transmission probability 
becomes zero. 



In contrast to the kicked rotor, the kicked Harper is 
found to exhibit a series of breakups and reformation of 
KAM barriers. These transitions, quantified by the expo- 




FIG. 2. The diffusive regime in 2-d parameter space . The 
unshaded part (top) in the is the regime of KAM barriers for 
momentum transport with /?ci ~ 0. The corresponding quan- 
tum plot(bottom) with rih — 300: In the unshaded part, the 
transport exponent is /3 < 10. 

Important feature of kicked systems with toroidal 
phase space are the AMs which are regular (stable) space- 
time structures coexisting with the chaotic sea in phase 
space and are accompanied by an hierarchy of island 
chains inducing anomalous transport /3ci > 1. Fig- 3 
shows one such super-diffusive parameter window whose 
origin is traced to a period-8 AM ||lj]. Once again, the 
quantum state lu = although localized exhibits a very 
strong enhancement of localization length in the classi- 
cally super-diffusive regime. It should be noted that in 
contrast to the diffusive peaks, super-diffusive spikes are 
in fact groups of many spikes exhibiting sensitive depen- 
dence on the parameters and hence describe transport in 
fractal phase space. 
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FIG. 3. For a fixed K = 0.86 the figure shows classical 
anomalous transport due to AM. In the corresponding quan- 
tum results, the transmission probability was zero except at 
the crosses for rih = 300 and the lattice size A'' — 43289. More 
precisely, the crosses indicate points where localization length 
varies between 400 to 800 which is about ten fold increase from 
the locahzation length in the diffusive regime. 

It is rather surprising that a pure quantum state uj — 
|12| can exhibit such a clear signature of almost all the 
classical transitions. This property may be associated 
with possible special structure of a; = state related to 
the fact that its quasi-energy is constant and thus insensi- 
tive to variations of parameters. We would like to point 
out that we have only investigated L > K part of the 
parameter space whereas duality implies that analogous 
behavior will be seen for K > L in x-space. 

In our earlier studies M for h w 0(1), we found patches 
of ballistic (localized) regions ior L > K {K > L) [|j. 
Numerical studies for smaller values of Ti suggest that 
the overall measure of the extended (localized) regimes 
for L > if [K > L) approaches zero as ?i ^ 0. 

We now discuss the symmetric Harper model with 
K = L. Here the quasi-energy states remain critical and 
hence exhibit diffusive transport for all values of the kick- 
ing parameter K . As ?i — *■ 0, (see Fig. 4) the transmis- 
sion exponent becomes fi independent provided we use 
renormalized parameter K instead of the bare K. The 
model exhibit transmission characteristics of the Harper 
equation for K < tt/2. Precisely dX K = 7r/2, the trans- 
port exponent begins to exhibit an oscillatory behavior 
(with frequency proportional to Ti). As opposed to clas- 
sical mechanics, where infinitesimal perturbation leads 
to chaotic regions whose size increases as the perturba- 
tion increases, the perturbation of such quantum system 
causes no immediate change in the transport character- 
istics. This suggests that roots of these transitions may 
be topological jlj]. It should be noted that the onset 



to oscillatory behavior is seen at higher odd multiples of 
7r/2, however, the behavior becomes prominent only at 
very small h. Other fascinating feature is that beyond 
K — TT, transmission probability appears to be periodic 
in K with period tt. Finally, the model exhibits a series of 
"transitions" precisely at K = Itt,1 = 1, 2... characterized 
by a discontinuity in the transport exponent. 
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FIG. 4. Variations in (3 characterizing the transmission 
probability along the line A" = L in kicked Harper model for 
rih = lOOOO(top), rih = 250(middle) and rih = 32(bottom). 
The plateau for K < 7r/2 is the transport exponent for the 
Harper equation. Although the exponent /3 oscillates with the 
N, all RG iterates show qualitatively the same behavior. 



Recently, a semiclassical analysis |11 hinted a possibil- 
ity of a series of enhancement of transport at ^ = Itt/2 
associated with the emergence of new periodic orbits of 
frequency I /A. The integer I is a kind of winding number 
that unfolds lu (confined to the interval 0— 2tt). The fact 
that K and not K determines the thresholds for various 
transitions poses a serious problem about the classical- 
quantum correspondence. As discussed in earlier studies 
pT| the symmetric model exhibits super-diffusive clas- 
sical transport at various critical values of the kicking 
parameter K. The most important of those are the 
period-1 and period-2 AMs which respectively occur at 
even and odd multiples of n. It is therefore tempting 
to associate K with the classical K as was suggested in 
Ref. ||l^. However, this singular scaling which associates 
the discontinuities in Fig. 4 with the classical AMs is in- 
consistent with the asymmetric case where no scaling of 
the parameters is needed to establish quantum-classical 
correspondence. It is possible that the transitions seen 
in the symmetric case are purely quantum-mechanical 
in nature and may have their origin in resonances of 
the (driven) RG flow and/or topological changes. We 
should mention that a possibility of some sudden changes 



at multiples of tt/2 also emerged in our analysis of the 
scattering problem. It turns out that the number of in- 
dependent propagating solutions of the scattering prob- 
lem [y, i.e. dimension of the S-matrix, increases by 1 
at Ki — Itt/2, I — 1,2..., matching with the discon- 
tinuities of the transport exponent P{K). Finally, the 
most challenging result which defies our intuition, is the 
(asymptotic) periodicity of the transmission probability 
2^lif+7r = T\i^, for K ^ TT which rules out any possibil- 
ity of quantum manifestations of classical super-diffusive 
transitions. 

Symmetric Kicked Harper model is an interesting ex- 
ample of a non-integrable system where the classical as 
well as the quantum transport is diffusive. It is in sharp 
contrast to the asymmetric case where the classically dif- 
fusive behavior corresponds to localized quantum trans- 
port. In view of this, it is rather surprising that in the 
asymmetric case the quantum system appears to respond 
to all the changes in the classical behavior, while in the 
symmetric case it remains insensitive to the variation in 
classical transport and instead repeats its behavior at ev- 
ery multiple of tt. This adds a new puzzle to the field of 
quantum chaos. 

Inability of the quantum system to delocalize in classi- 
cally diffusive regime and mimic the classical behavior for 
arbitrary small value of h as confirmed by RG analysis, 
remains an open frontier. Earlier studies have suggested 
phase randomization ^ due to classical chaos as a mech- 
anism for quantum dynamical localization. The fact that 
the kicked Harper model can exhibit localized, ballistic 
and diffusive transport irrespective of the fact that the 
corresponding classical system is chaotic challenges the 
phase randomization as the underlying mechanism for 
localization. Here we would like to propose an alterna- 
tive mechanism of dynamical localization: we speculate 
that the dynamical localization may be due to the can- 
tori barriers. These are invariant quasiperiodic trajecto- 
ries with infinite number of steps and provide an effective 
non-analytic quasiperiodic potential. The possibility of 
localization in quasiperiodic potential with infinite steps 
has been discussed recently [n5i. We would like to em- 
phasize that the scenario for localization due to cantori 
suggests that these barriers continue to inhibit transport 
even as ?i -^ irrespectively of the fiux through the holes 
in cantori. This scenario not only explains dynamical lo- 
calization in kicked rotor and Harper model (for L > K) , 
but also accounts for ballistic transport in kicked Harper 
for K > L. Furthermore, it is consistent with the dif- 
fusive quantum transport for K = L case , since the 
symmetric Kicked Harper model does not possess global 
cantori barriers (at least not of the type of broken KAM 
barriers) . We hope that further studies will put our spec- 
ulative views on solid footings. 
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